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SUMMARY 



Pulsed laser polymerization (PLP) is widely employed to measure propagation rate coefficients kp in free 

radical polymerization. Various properties of PLP have been established in previous works, mainly using 
numerical methods. The objective of this paper is to obtain analytical results. We obtain the most general 
analytical solution for the dead chain molecular weight distribution (MWD) under low conversion conditions 
which has been hitherto obtained. Simultaneous disproportionation and combination termination processes 
are treated. The hallmarks of PLP are the dead MWD discontinuities located at integer multiples of 
no = kptoCu, where to is the laser period and Cm is the monomer concentration. We show that chain transfer 
reduces their amplitude by factors e~'^*^^"" , consistent with numerical results obtained by other workers. 
Here Ctr is the chain transfer coefficient and Luq {L =integer) are the discontinuity locations. Additionally, 
transfer generates a small amplitude continuous contribution to the MWD. These results generalize earlier 
analytical results which were obtained for the case of disproportionation only. We also considered 2 classes 
of broadening; (i) Poisson broadening of growing living chains and (ii) intrinsic broadening by the MWD 
measuring equipment (typically gel permeation chromatography, GPC). Broadening smoothes the MWD 
discontinuities. Under typical PLP experimental conditions, the associated inflection points are very close 
to the discontinuities of the unbroadened MWD. Previous numerical works have indicated that the optimal 
procedure is to use the inflection point to infer kp. We prove that this is a correct procedure provided the 

GPC resolution a is better than nj^^. Otherwise this underestimates Lno by an amount of order a'^/no- 



Keywords : molecular weight distribution / propagation rate coefficients / pulsed laser polymerization / radical 
polymerization / theory 
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1 Introduction 



A fundamental material parameter in free radical polymerization (FRP) is the propagation rate constant^EEl 
kp governing the sequential addition of monomers to the active free-radical ends of "living" chains (see fig. ^ . 
In this paper we consider theoretically "pulsed laser polymerization" (PLP), ^ possibly the most accurate 
method to to measure kp. In PLP new radicals (living chains of one monomer, N — 1) are generated 
by a sequence of laser pulses which photocleave photoinitiators in the reaction mix. The laser is flashed 
periodically (the duration of each flash is extremely short (see fig. Ela))) with flash period to ~ O.lsec 
optimally chosen to be somewhat shorter than the average living chain lifetime. Thus every to the living 
chain population is updated by the injection of new primary radicals. After many cycles a periodic state 
is established where the time averaged rate of radical production balances the time averaged rate of living 
chain termination: pairs of living chains terminate to generate either a single ("combination") or a pair 
( "disproportionation" ) of "dead" chains, the final polymer product (see fig. E}. 

Now unlike the dead molecular weight distribution (MWD) of standard steady state FRP which has no 
singular featu res, in the case of PLP the dead MWD, (t)dead{N), is predicted theoretically to possess finite 
discontinuities'^'^ at chain lengths N = Lno, where L = 0, 1, 2, ... and 

no = Vpto , Vp = kpCu ■ (1) 

Here no is the number of monomers added to a growing living chain in time to and the "propagation velocity" 
Vp, namely the rate of monomer addition to a growing chain, is proportional to the monomer concentration 
Cm- The origin of the dead MWD discontinuities is the coherent sudden increase in termination rates due 
to the sudden increase in radical population following each pulse. In principle the discontinuities can be 
seen experimentally after measuring the resulting MWD by gel permeation chromatography (GPC). Thus 
no and hence Vp (or equi valently k^) i s inferred. This method has been employed to deduce kp for a variety 
of polymerizing systems. I ^ * ^ * ^ * ^*^ * ^^ * ^^ * 

In reality there exist a number of practical difficulties. (1) Broadening effects. Consider a living chain 
created at f = 0. After time to its length is not exactly Vpto] since polymerization is a random process, Poisson 
fluctuations arise around this length. This then leads to a softening of the theoretically infinitely sharp dead 
MWD discontinuities. Another, but quite different effect, is the direct broadening of the dead MWD due 
to intrinsic uncertainties in the measuring method. GPC equipment if supplied with a discontinuous dead 
MWD will output a reading indicative of a somewhat broadened discontinuity due to fluctuations in the 
transit time of a chain of a given length through the gel. (2) Chain transfer^ of radicals from living chain 
ends to surrounding unpolymerized monomer (fig. quantified by the parameter A. This is simply related 
to the chain transfer coefficient Ctr- 

ktrCx A 

= -rTT- = — . (2) 

"'p'-'M 'Jp 

where ktr is the 2nd order rate constant governing chain transfer to species X whose concentration is Cx 
(X is usually the monomer species). One expects chain transfer, important for species such as vinyl acetate 
and polystyrene, ^^21 to dephase the living MWD which would otherwise be a sequence of in-phase pulses. 
This will then affect the discontinuities in the dead MWD. 

Typical experimental MWDs produced by PLP are shown in fig. ^ where an obviously broadened first 
discontinuity is evident near molecular weight 40,000. In order to deduce the correct propagation rate 
constant, it is essential to have an algorithm to infer no from such a broadened discontinuity. Taking no 
as the molecular weight corresponding to the lower point of inflection, near the middle of the broadened 
discontinuity, has been proposed as the optimal approach. ^ Other authors have advocated that under 
particular conditions this criterion should be replaced with the maximum point. i^^*!^ 

In order to design and interpret PLP experiments many researchers have employed numerical simulations 
of PLP chemical kinetics. By varying the parameters of the simulations the resulting dead MWD is modeled 
under different experimental conditions. In such numerical studies Poisson-broadening is automatically taken 
into account and numerous side-reactions such as chain transfer can be incorporated. Guidelines can then 
be derived for extracting kp from the broadened MWD peak. Such numerical approaches have so far been 
employed both in the absence^^EH and the prese nce of chain transfer. DSEIESI xhe effects of experimental 
broadening have been addressed numerically in refs.^ESl 

Analytical solutions of PLP dynamics are also very desirable. Where available they provide a systematic 
and quantitative framework which exactly articulates the conclusions of the underlying physical model. 
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beyond the scope of numerical analyses. Analytical expressions can also greatly facilitate experimental design 
by obviating the need for numerical work. The first analytical solution of PLP dynamics appeared in refs.'^'i^ 
and included the effects of Poisson broadening for termination by both combination and disproportionation. 
A new derivation and a review of these theoretical results is given in ref.^U A closed form for the dead 
MWD for arbitrary sequences of laser pulses was derived in r ef.*^ Analytical work including chain transfer 
exists for the case of termination by disproportionation only. 122121 

In this paper we develop a completely general analytical theory for PLP, including the effects of broaden- 
ing and chain transfer analytically. We thus quantify the effect of these two mechanisms on the theoretically 
infinitely sharp dead MWD which enables us to discuss the validity of empirical rules on deriving kp from the 
broadened MWD peak. We consider a general (and typical) situation where both combination and dispro- 
portionation events occur. This extends previous analytical studies. We derive the most general analytical 
expressions which have been obtained hitherto for PLP, to the best of our knowledge. Closed expressions 
are derived in terms of the following independent parameters: Vp, kc, kd, A, to, and the concentration of 
radicals produced per laser pulse. Here kc and kd are the combination and disproportionation rate constants, 
respectively (see fig. EJ- An important quantity in what follows is the net termination rate constant, 

kEE kc + kd ■ (3) 



Our analysis describes low conversion FRP which allows us to make the approximation that the net living- 
living bimolecular termination rate constant, k, is independent of chain length (see fig. [SJ. In reality k{M, N) 
depends on the degrees of polymerization M, N of the reacting chains. However in dilute solutions (under 
good solvent conditions) theory ^'"^ predicts a very weak dependence on chain lengths: k{N, N) ~ N^"" 

with a ~ 0.16. A very weak dependence of k on chain lengths has also been well established experimentally 
in phosphorescence quenching studies. ^^Sl At low conversions, living chain terminations occur in a pure 
monomer solvent (usually a good one) and thus the approximation of constant k is expected to yield rather 
accurate conclusions. Making this approximation, all quantities can be calculated exactl y. The effect of the 
dependence of termination rates on chain length has been addressed numerically in refs . | 29 | 30 | 3i | 32 | 33 | 34 | i8 | 
which employed various empirical rules for the form of k{M,N). 

The living MWD, (p{N,t), like all quantities settles down to a periodic form with period to. Bimolecular 
reactions as in figs. |21 El generate dead chains of length at a rate 

I 

<^dead{N,t) = -kc / dMcl){M,t)cj,{N - M,t) + kd<j){N,t)^{t) + Xcj>{N,t) , (4) 
^ Jo 

where ^ is the total concentration of living chains. The resulting dead MWD will include Poisson broadening 
effects provided one uses the correct Poisson broadened living MWD in eq. Q). In the following we will 
employ the notation 4>dead to denote the dead MWD without Poisson broadening, including however chain 
transfer effects (i. e. (pdead results from using in eq. Q the form for 0(Af, t) which correctly includes chain 
transfer but which neglects Poisson broadening). The full Poisson broadened MWD is denoted 4'^J°^''f while 
the reading output by the GPC measuring equipment we name (jidtad- 

The dead MWD is calculated in successive steps as follows. In section 2 we review "ideal" PLP, i. e. we 
ignore chain transfer (setting A = 0) and broadening effects. The time averag ed dead MWD, turns 
out to have finite discontinuities at A^ = jno {j is a positive integer). In section 3 we generalize the results of 
section 2 to include chain transfer and express the resulting dead MWD, (pdead, in terms of (f>^ffi^d- We find 
that chain transfer diminishes the amplitude of the dead MWD discontinuities, but does not broaden them. 
Broadening effects due to both Poisson broadening and due to limited experimental resolution are included 
in section 4. We show that the inflection point rule is correct for inferring hq only when the width of the 
Poisson broadened discontinuity is larger than the experimental resolution. We conclude with a summary 
of our results in section 5. 

Throughout, [N] will denote the nearest integer to N/uq which is less than N/uq, whilst rN = rem[A^, uq] 
denotes the remainder of N modulo uq. Similar notation will apply to time i, with uq replaced by to. Thus 

N=[N]no + rN, rN = Tem[N,no]; t = [t]to + rt , rt = rem[t,to] . (5) 



2 Ideal PLP 



Let us first briefly revi ew th e "ideal" PLP problem ignoring the effects of broadening and chain transfer, as 
analyzed by Olaj et alJ ^ ^ ^^I Termination may occur either by combination or disproportionation reactions. 
First we calculate the living MWD which depends on the total termination rate constant, k, only, and not on 
the mode of termination. We then derive separate expressions for the dead chains terminated by combination 
and disproportionation, respectively. 

2.1 Total Concentration of Living Chains ^(t) and Reaction Field H{t) 
The number density MWD of living chains, (j>{N,t), obeys 

</) = ~vpq^' ~ B{t)<t> + 6{N)R,{t) , (6) 

with (\){N < 0,i) = 0. Here = d^jdt and 0' = d^jdN. Note that the S{N)Ri{t) term which injects 
primary radicals (chains of length = in our continuous framework) into the system at rate Ri(t) per unit 
volume, is equivalent to the boundary condition Vp(f>{0,t) — Ri(t). In PLP dynamics the time-dependent 
primary radical production rate, Ri{t), is a train of pulses, 

oo 

R,{t) = R,to '5(^+.?^o) , (7) 

j=~oo 

as shown in fig. |SJa). Ri is the time averaged radical production rate. The —Vpcj)' term in eq. ^ describes 
propagation whilst H{t)(p{N,t) is the total termination rate (due to both combination and disproportiona- 
tion) for a chain of length N due to the "reaction field" H{t), namely the total termination rate per unit 
time for a given chain due to all other chains in the system: 

/•oo 

H{t) = k'i>{t) , *(i)= / dN(l){N,t) . (8) 



Here ^(i) is the total number density of living chains whose dynamics follow after integrating eq. © over 
all positive N values: 

= -Hit)-^ + R,{t) ^ -k'^^ + R,{t). (9) 

For t different from integer multiples of to one has Ri (t) — and thus 'J obeys simple second order reaction 
kinetics, ^I^ — — fc^I'^, with solution 

*W = 1 (0 < t < to) , (10) 

1 + <KWo 

where ^'o = ^(0+) is the value of ^' just after a new laser pulse at t = 0. We define (3 as the surviving 
fraction just before the next pulse at to-, 

P^^{to)/^o . (11) 

The values of "ifo and /3 when the system settles down to a periodic stationary state are determined in terms 
of Ri,k, and to, by equating the number of terminated chains in one cycle with the number of new living 
chains introduced by a pulse, 

(1 - /3)*o = i?<io ■ (12) 
This relation together with eqs. (|10|l and Hll|) yields 



1 - 

*o = ^Rih 



l + {l + 4/(i?,H2)}i/2l , ^ . (13) 

J 1 + fcWoto 



The steady state vE'(t), sketched in fig. Efc), is of course periodic with period to, and has discontinuous 
jumps AvP = (1 — /3)^o at integer multiples of to- 



A 



2.2 Living MWD 



The solution to eq. © is derived in appendix A for a general Ri{t) in terms of H{t) whose important 
properties are derived in appendix B. Using the present periodic form for Ri{t), eq. using eqs. (|12|) and 
(|B2|I and recalling the periodicity of \1/ and H, we have 

oo 

= <l>env{N)S{N - Vpt ~ jno) , (14) 

j = -oo 

where 8 is the step function (8(x > 0) = 1, 8(.t < 0) = 0) and the function (j)env, sketched in fig. Efb), is 
the living MWD envelope 

</.e„„(iV) ^ 8(A^)/3[^l(l-/3)vl/(7V/^;,) = e(A^)/3[^l(l - /3)^^^|i;^ . (15) 

The living MWD, sketched in fig. Elb), is a train of pulses in phase with Ri{t). A new 5-pulse is born 
at = at time t = 0, say, and then moves in the +A^ direction with velocity Vp. When in the interval 
(j — l)rto < N < jriQ (the "jth sector") this pulse suffers the same field H{t) as in all other sectors. The 
amplitude of the i5-pulse in sector j + 1 is thus reduced by a factor /? with respect to sector j and follows the 
envelope (t>env{N). 

2.3 Dead MWD 

—rideal 

The time-averaged rate of dead chain generation, 4>ciead i-^) i the sum of the contributions due to termination 

-rideal —comb -rdisp 

by combination and disproportionation, ip^f.ad (-^) — 4>dead i-^) + 4'deadi^)- Each of these two contributions 
can be found by inserting eq. (|14|l into the first and second terms on the rhs of eq. I@J, respectively, and 
then time averaging over one period. In the following we will use however a simpler method. 

Consider first chains terminated by combination. We notice that, due to periodicity, the total concen- 

—comb 

tration of dead chains of length N produced by combination per period, (j^dead (-^)^0i equals one half the 
concentration of dead chains produced by combination by one single S pulse throughout its entire lifetime. 
Thus choosing the j = pulse in eq. H14|l . bimolecular kinetics imply 

—comb k r°° 

(^dead (^) = ^ / '^^^ 4>env{M)S{M - Vpt)(f,{N - M, t) . (16) 



J- 



oo ^ — oo 



This is proved rigorously in Appendix C. Since the integrand of eq. (|16|l is zero for negative t and M values, 
we have extended the lower limits of the integrals from to — oo for future convenience. 

Substituting the expression for (f> of eq. (|14() into the above expression, it is shown in Appendix D this 
leads to 



where 



^dead (N) = Cf\{L + 1)/?-/^ ^ ) + Lf3--'f' ( ) } , (17) 

N^Lno + rN, /(n) ^ \ , ^ ^ ^^fi^ . (ig) 

fj + [I - (j)n/no Avpto 



Similarly, the time-averaged concentration of dead chains generated by disproportionation can be cal- 

-rdisp 

culated by noting that 4'deadi-^)^o is equal to the concentration of dead chains of length N produced by 
disproportionation by one single (S-pulse throughout its lifetime. Choosing j = for this pulse in eq. (|14|l 
one has 

—disp u, r°° 

<f>deadiN) = / dt <t>enAN) <5(iV - Vpt) ^ {t) . (19) 



Integrating the 5 function in eq. (|19|) and using eq. H15|) to express (penv in terms of 5* yields 

0rfead(^) = J±-f(rN) , N^Lno + rN . (20) 



k {l-(3' 

-rideal —rcomb —rdisp 

The dead MWD, = tc/)^^^^ = t{(t>^^^^ + (jiaead)^ is shown in fig. Hd). Both combination and 

disproportionation contributions have discontinuities at the end of each sector L (where L — Q labels the 
first sector, < < tiq)- The discontinuity, ^4>^^^ad^ boundary of sectors L ~ 1 and L is 



3 Chain Transfer 

In this section we will study how chain transfer modifies the results of section 2. Each transfer event results 
in one dead chain of length N and one living chain of length TV = as depicted in fig. O 

3.1 Living MWD 

The living dynamics now include additional loss through transfer from chains of length N . These radicals 
are reinjected as a source term at iV = 0: 

0=-t,p0'-[H(t) + A]0+[i?,(t) + A^'(t)],5(iV) . (22) 

The total reinjection rate is A^* where A is defined in terms of the chain transfer coefficient in eq. 1^)1. 

Integrating eq. H22|l over all positive N values (0 is zero for negative iV) one immediately sees that the 
total number density of living chains ^'(i) obeys the same dynamics as without chain transfer, eq. ©. This 
must be true, of course, since transfer events conserve the total number of living chains and k is independent 
of chain length in the present approximation. The crucial point is that ^{t) and Hit) = —k^(i) appearing 
in eq. (|22|l are the same periodic functions of time as without chain transfer, eq. I|in|l . 

Since Hit) is a known function of time, the dynamics of eq. 1)22(1 are effectively linear and we can use 
the superposition principle to write = 4>coh + 0mc where 

4>coh - -M'coh - + A<^coh + Mt)S{N) 

^^nc = -Vp^[„^-[H{t)+\]4>rac + >^^{t)5{N) . (23) 

Here 4)coh is the coherent part, which as we will immediately see is a sum of pulses in phase with radical 
production, whilst (j)inc is the incoherent part which is out of phase due to random transfer events. 

Replacing — > + A in the general solution, eq. I|A3|I . to the living dynamics without chain transfer, 
one sees 4>coh is e~^^/'"p times its no-transfer form of eq. (|14|l . Similarly replacing H H + \ and Ri —> A^, 
one obtains (pine after using eq. I(B2(I . One finds: 

°° (l_ g)e->'t 

Cpcoh{N,t) = e-^^/'"" ^ 0e„„(A^)5(A^- V-i"0) , ^coh(0 - \ _ ^j_Ato ' 

j=-oo 

\ -XN/vp 

<t>^nc{N,t) = Vl/(t)/3W~[*~^'/''ple(iV) , ^^nc{t)^^{t)-^coh{t) , (24) 



P 



where \E'coh, ^inc are the corresponding total concentrations of chains in each part of the MWD. (pine is a 
series of propagating pulses each of width no, the amplitude of which decreases with molecular weight and 
time. Notice however that unlike (pcoh, <t>inc extends over all N values at a given time, (pine has discontinuities 
at iV = Vpt + jriQ. At each discontinuity, (pine decreases suddenly by a factor /3. 



3.2 Dead MWD 



Dead chains are generated by both Uving chain termination (first two terms on the rhs of eq. 10}) and chain 
transfer events (last term of eq. 10} ). We consider first the contribution due to termination by combination. 
Using the results of Appendix C, the time averaged rate of production of dead chains by combination is 
given by 

—comb k f°° f°° .,r, 

^dead {N}^^ dt dMe-^^''^-4>env{M)5{M - Vpt) {q^,oh{N - M, t) + 20„,(iV - M, t)} (25) 



2to J- 

This is the contribution to (t>dead ^o termination by combination by one single d pulse belonging to the 
coherent living MWD, throughout its lifetime. Combination occurs either with other coherent chains or by 
"cross-coupling" with chains belonging to 4'inc- In eq. H25() we have assumed that chain transfer is a weak 
effect, i. e. that A is small in a sense to be quantified below, and we keep the leading correction only. Thus 
we discarded the 0{X^) 4'inc4'inc termination term. 

Now the coherent-coherent coupling term in eq. H25|l is equal to that of eq. H16|) . multiplied by q~^^/^p _ 

—rcomb 

The only difficulty is the coherent-incoherent cross-term (^^g^j^; cross- Appendix E it is shown that 
—comb Ae^-^^/^p k ( r r /rN 

'l^dead.crossiN) = " ^^^^ + " ^l^" + + " 

Vpto k I \2vp 



L(3'^~\l-(3)m\——^j-(3^^>\—j^ . (26) 

—disp —trans 

Let us calculate now the disproportionation and chain transfer contributions, (t>dead ^dead ■ Time 
averaging the last two terms of eq. (0J one has (bar denotes time average) 

—rdisp —rtrans 



'/'deadW+'^dead (^)=fc<i0(^)* + A0(iV) ■ (27) 



Now from eq. (|22|l one has k(j)[N)'^ + \cf){N) = —Vpd<f)/dN, using (f> — Q in stationary state. Thus using 
eq. (123, 

<l>dead{N) + C^dead {N) = "^^P^ + f A0(iV) . (28) 

The rhs of eq. (PS} is calculated in Appendix F (see eq. HF5|) ') using the expressions for (j) derived in section 
3.1. 

— —comb —disp —trans ^ 

Thus collecting all terms, 4>dead — ^dead + 4>dead + 4'dead i ^^d noticing that the kc/k term of eq. I|F5|I 
cancels the last term of eq. (|26|l , we have 

— , (—ideal -^inc "] A\ AXBfn 

^dcadiN) = e-^""/^^ [N) + e </.,,,,(A^) j , e ^ (J^^y^ = ' ^^^^ 



where 



^dcadiN) ^ C 



k ( / rM 

^ |/3^[(L + 1)13 -L]+I3\L + 1)(1 - 0}f ^— 



(30) 



The two contributions to the dead MWD, (pdead — t4>deadj ^'^e plotted in fig. |^ The incoherent part due to 

— mc —inc __ 

transfer events has no discontinuities, i.e. 4'deadi-^^o) — <^(iead(-^^o ); may be seen using eq. H18(l . It is 
proportional to the small parameter e ~ XIq. 

—^deal 

■~> ' : ■ : the only discontinuous term in eq. (|29|l is 4>dead j using eq. (|21|l the discontinuity at iV = Lng is 

AcPdeadjLno) ^ ALto nL 2 + ip 

0dead(O) l + 4p/(l-/3)+e(l + 2p) ' ^ 

where A(l)dcad{Lno) = (j)dcad{Ln^) - 4'dead{Ln'^). 
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4 Broadening and the Inflection Point Rule 



The MWD derived in the previous section, eq. (|29|l . has ignored (i) Poisson broadening and (ii) broadening 
due to the hmited accuracy of GPC measurements. 

Consider first Poisson broadening. The true length distribution of a chain growing with mean velocity 
Vp is not 6{N — Vpt) as has been assumed in sections 2 and 3, but rather the Poisson distribution (see fig. d) 
which for large N tends to a Gaussian: 

P^(Ar) ^ l!;^e-"-* « ce-(^^"''*)'/(2-p*) ^ (TV » 1) . (32) 

Here c = {2'KVpt)~^/^ is a normalization constant. 

Now since the termination rate constant k is independent of chain length, clearly broadening leaves the 
total number of living chains ^'(i) unchanged. Hence the field H{t) is unchanged, and thus a group of living 

~ f* H . 

chains injected at t = are still depleted in number by the factor e Jo after time t, regardless of which 
chains are a little longer than Vpt and which a little shorter. 

Moreover we show in Appendix G that even in the presence of chain transfer, the Poisson broadened 
living MWD, (f>broad{N), can be expressed in terms of the unbroadened living MWD, <j>{N), calculated in 
the previous sections as follows: 

/OQ -1 _ 

dMcj^iM, t)A{N; M) , A{N; N) = ^^^-(^-Nf /(2N) (33) 

-00 V27r7V 

Thus the (5-pulses of in eq. H14|) (or of (pcoh in eq. I|24|) if chain transfer is present) are replaced by 
A's, namely Gaussians of mean N and width VtV. Note that with increasing N , the relative width of the 
broadened A-functions decreases as I/VtV. 

Then substituting eq. H33|) into eq. Q we find that the dead MWD, (j>^J^, is also obtained by integrating 
the dead MWD in the absence of Poisson broadening, eq. (|29|l . against a A function: 

/oo 
dM^dead{M)A{N;M) . (34) 
-00 

(We used the identity J^^ dMA{M;P)A{N - M;L) = A{N; L + P) for the termination by combination 

term of eq. Ql). Here 4)dead = t4>dead- example of the broadened MWD is shown in fig. |S1 where 
't'^dead Calculated by numerical integration of eq. H34|l . Eq. H34|) has been called "a posteriori Poisson 
broadening" In the present work we proved its validity. 

Now the accuracy with which the resulting dead MWD of eq. (|34|l is measured is limited by the resolution 
of the GPC equipment. Thus the MWD implied by the GPC output reading is 

/>oo 

read / T\T\ / A^^oad 



fdTadiN) = / dM c^TaTm G{N- M) , (35) 

J —00 

where G{N; M) is the experimental reading (normalized to unity) given a perfectly monodisperse MWD 
input of chains of length M. 

Consider values of N near the discontinuity of the unbroadened dead MWD, N = Luq, and let us examine 
the location of the inflection point. In the following we will assume that near the discontinuity to leading 
order G can be approximated by a Gaussian distribution, G{N;n) = {2Tra'^)~^^'^e~''^~"'^ /^^'^ \ where a is 
a measure of the experimental resolution near Lno. Substituting <j>''Jf^f from eq. JSl in eq. (|35|l one has 

/.oo _ _ -iN-Mf/l2{M+a^)] 
4>7ead{N) = / dM <j,aead{M)A{N- M) , A{N , M) EE — — — , (36) 

after using dP A{P; M)G{N; P) = A{N; M) for Gaussian G. 

In Appendix H the curvature of 0^g°5^ is evaluated close to a discontinuity (eq. (|H4() ) . Setting this to 
zero one finds two roots, one irrelevant at 5N = N — Luq « 2{Lno + cr^), and one corresponding to the point 
of infiection at 

6N « j^+-^(Lno + , 0+/- ^ ^dt1iL4^-) , ^'+/- ^ ^^{Lnf-) , (37) 



Now in PLP experiments, parameters are normally chosen such that the dead MWD decreases by a factor of 
order unity between successive peaks. Hence cf)'^ is of order (j)^/no and similarly for Since (/)'(_ < 4>'_ < 
and > (j)- > 0, and taking L = 0(1), one has from eq. (|S7|l . 

SN « -A{1 + Ba^/no) , (38) 

— 1/2 

where A,B>0 are constants of order unity. This theoretical result is valid provided Uq <C 1 and /? is of 
order unity. Note that SN is always negative. 

1/2 

Consider first very good experimental resolution, a ^ Uq . In this case the observed broadening is 
predominantly Poissonian. According to eq. I|38|l . the displacement of the inflection point of the measured 
MWD away from Lhq is then of order one monomer. This is an unimportant correction which is of the same 
order as the accuracy of the continuum description of the MWD employed in our analysis. We conclude that 
for this case the inflection point is the appropriate point from which to extract kp. 

1/2 

Now when a ^ , the broadening of the measured MWD is dominated by the finite experimental 
resolution. In this case the inflection point is of order a"^ /no 3> 1 monomers away from Ltiq. Therefore 
for large a, using the location of the inflection point to deduce no leads to a systematic underestimation 
of Vp. We remark however that this is more accurate than using the local minimum or maximum of the 
MWD which will then be at a distance of order cr away from no (since we have assumed that cr <C noi hence 
cr > cr^/no). 

As a specific example, consider the experimental MWD of fig. 2| In this figure, the distance between 
the local minimum and maximum of the MWD around the inflection point near molecular weight 50,000 
(no w 500) is of order 200 monomers. In this case Poisson broadening would only account for a width 

1/2 

« n-Q « 22 monomers. Hence broadening appears to be mainly due to limited experimental accuracy and 
cr sa 100. The relative error using the inflection point for the calculation of Vp would then lead to an error of 
order cr^/nQ w 4%. 



5 Conclusions 



In this work, we have calculated analytically the number MWD of dead chains produced by pulsed laser 
polymerization at low conversions. The following effects have been included: living chain termination by both 
combination and disproportionation; chain transfer of living chain radicals to the environment; broadening 
effects. Ignoring broadening, we have found the MWD is given by 



4A/3io 
(1^ 



(39) 



where 
'AS(A^)= const. 

^deUN) ^ const. 



/(") 



k 



-j\rN) 



/3^[(L + l)/3 -L]+ p'^iL + 1)(1 - (3)f 



k 

+L/3^-i(l-/3)/ 
/3 + (1 - /?)n/no 



rN 



(40) 



Here N = Luq + rN, and "const. " is a normalization constant. The dead MWD is discontinuous at integer 
multiples of no and is a function of the following independent parameters: no, A^o, fcrf/fc, kc/k, and (3. Here 
P (see eq. (|13|) '). the fraction of living chains surviving between 2 successive pulses, is a function of the 
dimensionless parameter kRit^, where Rit^ is the concentration of new radicals produced per pulse. Eq. 
(|39|l applies in the limit in whi ch cha in transfer is weak, i.e. Ato <C 1. Setting k^ — 0, eq. H39|) redu ces to 
the dead MWD derived in refs. ,1^2131 while setting A = one recovers the dead MWD derived in refs.'^'^ 
An important quantity in eqs. H39|) and H40I) is the relative magnitude of the dead MWD discontinuities, 
shown in eq. 131() . Optimally large discontinuities typically correspond to values of /3 intermediate between 



zero and unity as has been thoroughly examined by numerical simulations^^^ and discussed in ref.^ For 
/3 <C 1 these relative magnitudes are undesirably small. The limit /? ^ 1 (realised when kRitQ is small) may 
be undesirable also, but this depends on the mode of termination: (i) For the extreme of termination by 
disproportionation only (p — s- cx) in eq. the discontinuities vanish in this limit while (ii) for termination 

by combination only (p — 0) the discontinuities attain their maximum value as /3 — > 1. Thus if combination 
is strongly dominant, by choosing small Ri this may be a desirable limit (though noise effects are then 
particularly strong since the absolute amplitude of the dead MWD is very small). In a typical case though 
where termination occurs by both combination and disproportionation the optimal choice corresponds to 
intermediate f3 values. 

The effect of chain transfer is to reduce the amplitude of the entire MWD. The amplitude of the dis- 
continuities is reduced by a factor Q-^toN/ng^ (j^ follows that if A^o > 1, most living chains suffer a chain 
transfer event before growing to a length of order Hq; consequently the MWD amplitude at TV = rig is small. 
Under such conditions PLP is not a useful method for measuring uq.) In addition to these effects, transfer 
introduces an extra "incoherent" contribution to the dead MWD (see fig. ^ which has no discontinuities 
and whose magnitude is a fraction Xto of the dead MWD. 

We have also treated 2 types of broadening: Poisson and distortion by MWD measuring equipment. 
Poisson broadening smoothes out the discontinuities of the dead MWD of eq. H39|l . The broadened dead 
MWD can be computed by performing numerically the integral of eq. H34|l . We found that for /? of order 
unity the position of the inflection point of the smooth MWD is not shifted significantly away from the initial 
discontinuity. ( For very small /3 the inflection point rule leads to underestimation of no.) Thus, if Poisson 
broadening dominates, then its location can be used reliably to infer kp. 

In practice, the 2nd type of broadening, due to limited GPC resolution, may swamp the Poisson effect. 
We found that when the band width a of the GPC output reading, given a perfectly monodisperse sample 

1/2 

of degree of polymerization Luq, is much greater than ng , experimental broadening then dominates near 
the discontinuity at Lhq. This leads to a displacement of the point of inflection by approximately cr^/no 
monomers below Lhq. We remark however that (for f3 of order unity and a <^ uq) the point of inflection is 
closer to Uq than the position of either the minimum or maximum peaks. 
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A Solution to Living Chain Dynamics for Arbitrary Ri{t) 

In this appendix we solve eq. © for general radical production rate Ri(t). Defining £^ = t — N/vp and 
/(^, t) EE (j){N, t) we have from eq. ® 

df 

-± ^ = -H{t)f + S{vpt ~ VpOR, it) , (Al) 
with (/)(7V < 0,<) = implying /(^ > t) = 0. The solution of eq. is 

fiU) = - r dt'e-I>d{t'-m{t') = e{t^O-e'^^"m) , (A2) 
where Q is the step function, Q{x > 0) = 1 and 9(x < 0) = 0. Hence 

<^(iV, t) = ^'-"/". R.{t - AT/,^) . (A3) 

Dp 



B Properties of the Reaction Field H{t) 

Since 'J — —H'^ within one cycle (see eq. (O), it follows that if both ti and t2 belong to the same cycle 
then 

g- 1: = vi/(i2)/vi/(ii) , ( [t,] = u ) . (Bi) 

Since ^'(t) drops by a factor (3 in one cycle of duration to we also have 

~ H{t')dt' ~ f^*^^'"' H{t')dt' ^ r''' H(t')dt' *(t2) r, 1 Ft 1 . , 



after using eq. (|B1(I . 

C General Expressions for Dead MWD's Generated by Combina- 
tion 

Consider two living MWDs (pi , 4)2 which are periodic in time with period such that (pi can be written as 

00 

01 (N, t)= Al {N)fi {N - vpt - jno) , (CI) 

j=-oo 

where fi has period uq. Suppose dead chains are generated only by bimolecular reactions involving one 1 
and one 2 chain. If only combination reactions occur, the time averaged rate of dead chain generation then 
satisfies 

'i^deadiN) = - dt dMMM,t)MN-M,t) 

^0 Jo J-00 

tQ poo 



- V / dt dMAi (M)/i (M - vpt - jno)(^2 (N - M, t) 

, 00 /■(i + l)to fOC 

- '/ dMAi{M)h{M - Vpt')MN ~ M, t' + jt^) 

''0 ■ „ Jitn J -00 



] = -oo 



(C2) 



1 1 



where we changed variables to Vpt' = Vpt + juQ. We remark that when there is only one chain type reacting 
with its own type to produce dead chains, then (f>i = (/)2 and the rhs of eq. ljC2p should be divided by 2. 
Since <f)2 has period tg one has the following equivalent expression: 
— u r°° r°° 

(l^deadW - - / dt' dM A,{M)h{M - Vpt')MN - M, t') . (C3) 



With Ai{N,t) (t)env{N), /i ^ (5, 02 ^ 01, fc — > fcc and dividing by 2 this yields eq. of the main text. 

D Ideal Case: Combination Part of Dead Chain MWD, eq. (17) 

Changing variables n = 2vpt + Juq in eq. H16|l and using eq. (|14|l to substitute for in eq. H16(l one has 

4>dead W = j dn ^ <j>env ^— j (f>env ^— j S{N - Tl) 



Avpto 



Writing N = Ltiq + rN, with L = [N], we evaluate the sum of eq. IjUll) by noting that since (p^nv is zero 
for negative values of its arguments, only —L < j < L contribute. Now using eq. (|15|l one has 

(L + j)no + WV\ ^ / {L- j)7io + rN \ ^ f l3''^ln.{rN/2) {L - j even) 

2 J 2 J [ ^L-i^2^^^((,Ar + ^„)/2) (i-jodd) 

{~L <j<L) 

(D2) 

Noting that the number of even and odd L~j values in the interval —L <j<L is L + 1 and L, respectively, 
and using eqs. (|D2|) and (|D1|I . eq. p7|) of the main text is derived after expressing 0e„„ in terms of 5* using 
eq. (dH). 

—comb 

E Derivation of Coherent-Incoherent CoupUng Term 0dead cross 

Using eqs. H24(l and eq. (|15|l . one has from eq. (|25|l after integrating over the (5- function: 

Vpto Jo ^ J 

= ^ (1-/3)/ dA/*2(M/z;„)/3^+i+['-^-2'-A^l 

Vpto Jo 

N = Lno + rN , M = Puo + rM . (El) 

Here we used the identity — [A^ — 2M] = [2M — A^] + 1. The integrand of eq. (|E1|) is independent of P and 
only depends on the relative magnitude of rN with respect to rAI. Since O < P < L we may split jj^ into 

L integrals /J^^^^ + /rjv/2^"°^^^ + /(r';v+rio)/2 corresponding to P < L — 1, plus the P = L term whose upper 
integration limit is rN instead of tiq: 

-comb Afcce"-^^/"" 

Vpto 

{L + l)/3^+i / d{rM)-i'^irM/vp) + {L + 1)(3^ / d{rM)^^{rM/vp) 

Jq JrN/2 

+L(3^-^ d{rM)-^^{rM/vp)+ d{rM)-^^{rM/vp) \ (E2) 

J(rN+na)/2 J{rN+n(,)/2 I 



t'dead,cross W = — ( 1 " /?) X 
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—comb __ ^ 

The final expression for <j)dead,crossj ^1- of tlie main text, follows from eq. HE2|) after using the identity 
~k f*^ = *(t2) - for < ti,t2 < to- 

F Disproportionation and Chain Transfer Contribution to Dead 
MWD 



From eq. H24|) . the time average of (f)coh is 



(b,ohiN)^^e-^''^^''^enviN) . (Fl) 
Vpto 

Hence using the definition of (j)env (eq. (|15(l ^ and = — fc^tf^^i) for < i < to, 

- %^ = ^e-^^/''''^,„„(7V) [A + fc*(r7V/«p)] (F2) 

Civ fpEo 

The time average of (j)inc is calculated after integrating the 3rd expression in eq. (|24ll from to to- The 
value of the integrand depends on the relative magnitudes of rN/vp and t: 

cb^nciN) = / dt^{t) + - dt^{t)\ . (F3) 

^^pto I Jo P JrN/vp I 



Hence neglecting terms proportional to A^, 



-i-Ae-^^/"-0e™(iV) + O(A2) 



(F4) 



Thus from eqs. (|F1|) - (|F4|) one has 

after using eqs. p5|) and p8|l and the definition of e in eq. (|29|) . Here we neglected ©(A^) terms. 

G Poisson Broadening of Living MWD: Derivation of eq. (33) 

In order to include broadening effects in the living dynamics, a "diffusion" term must be added in eqs. © 
and as follows: 

4>broa, = -V,^^ + - H{t)<f> + R.mN) , (Gl) 

where H{t) = H{t) + A and Ri = Ri{t) + A^(t). In the absence of the last two sink and source terms, 
eq. (IG1|I describes the evolution of the probability distribution function a unidirectional walk with step size 
unity and velocity Vp, in the continuum limit. The diffusion coefhcient for such a walk'^*^ is Vp/2. 

Now in the absence of the source and sink terms, the propagator of eq. IjGlj) for a pulse generated at 
= at t = 0, is simply A{N;Vpt). (Here we neglected the condition (j)broad{N < 0;t) = which implies 
that the propagator be zero for negative N values. However for Vpt > 1, the amplitude of A(7V; Vpt) becomes 
exponentially small for negative N. This implies that the use of A is a very good approximation provided 
we examine the living MWD for iV ^ 1.) Hence the solution of eq. (|G1|I is 

4>broad{N,t)^ I dt'e'^^'"R^{t')A{N;vp{t^t')) . (G2) 



This solution is identical to the solution of the unbroadened living MWD in eq. ljA2p . after replacing 
5{t' — = Vp6{N — Vp{t — t')) in eq. l|A2p by A(A^; Vp{t — t')). Hence eq. (|33|) follows immediately as one 
may check after replacing (p from eq. IIA2p into eq. (|33() and checking that eq. (|G2p is recovered. 



H Curvature of Broadened Dead MWD Near the Point of Inflec- 
tion 



Considering molecular weights near Luq, expanding (pdeadin) one has (see eq. (|37f) for notation) 
(/>(iea<i(n-) = (0+ + 0+(5n + ...) 6((5n) + (0- + (?!)'_(5n + ...) 6(-(5n) , 5n = n - Lno . 



(HI) 



Since the amplitude of the unbroadened dead MWD decreases smoothly by a factor of order unity between 
successive peaks, the magnitude of cf)'^ is of order 0+/no, and similarly for </>'_. Now since the Sn which 
will survive in the integration of eq. 134(1 around N « Luq is of order E = (Luq + a^Y^^, eq. (|H1I) is an 
expansion up to order S/ng. Approximating A with a Gaussian of fixed width, and keeping terms up to 
order 1/E (note E/ng > l/S): 



A{N;n) 



SN 



(27rS2)i/2 
N - Lno ■ 



-{SN-Snf/{2T,^) 



Sn Sn{SN ~ Snf 
2E2 + 



(H2) 



Substituting eqs. l|Hip and l|H2p in eq. H34|) we obtain the following approximation to the MWD output 
reading near Ln^: 



^dead 



(N) ^ 



(2.S2)i/2 J, d{Sn)i^c^^+Sn^<l>'^- — 



(/)+ (f>+{5N-Snf 



2S4 



-(SN-Sny/{2-S^) 



(27 

A straightforward differentiation of eq. I(H3(I gives 

1 



l2±read 
" Vdead 



(27rS2)i/2 



2S2 



1 + 2SN - 



{SNf 



-(<5JV)V(2S2) 



(H3) 



(H4) 
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Figures 



FIG. 1. Schematic of monomer addition to the free radical end of a living chain. 

FIG. 2. Termination reaction of 2 living chains of length N and M. Combination results into one dead chain 
of length N + AI while disproportionation results into 2 dead chains of length N and AI. 

FIG. 3. Schematic of a chain transfer reaction with monomer as the transfer agent. 

FIG. 4. MWD of poly(methyl methacrylate) produced at 40"G by PLP with to = .Isec as measured in lei.^ 
Horizontal axis is molecular weight. Arbitrary units on vertical axis. 

FIG. 5. (a) Rate of radical production as a function of time. Ri{t) is non-zero only during the short duration 
of the laser pulse the period of which is to. (b) Snapshot of the living chain MWD (piN, t) during < i < to 
in the ideal case (absence of chain transfer and Poisson broadening). Here 0o = 0(0). consists of a 
series of (5- function pulses whose amplitude follows the envelope (penv {N) . (c) Number density of living 
chains as a function of time (see eq. H1Q(I ^. (d) Dead MWD in the ideal case. <j>^£^^j is discontinuous at 
integer multiples of no. 

FIG. 6. Contributions to the dead MWD in the presence of chain transfer. The upper curve is the coherent 
contribution. The lower curve is the incoherent contribution, proportional to the small parameter e. The 
slope of the incoherent MWD is discontinuous at integer multiples of no- Values of parameters used: 
Vp = 5000s-\ kc = kd = 2 lO^Af-ifi-i, A = .5s-\ to = .Is, R^to = 5 IQ-^M 

FIG. 7. Snapshot of living MWD including the effects of Poisson broadening: the (5-bmctions of the living 
MWD of fig. Elb) are replaced by Gaussians. 

FIG. 8. Illustration of dead MWD before and after broadening. The solid line is the unbroadened MWD. 
The dashed line is the Poisson broadened MWD. The dotted line is the reading on the GPC equipment, 
given the dashed line as the input MWD. We assumed uniform resolution cr^ = 5000 (see main text). 
Values of parameters used: Vp = 5000s"^, kc = kd = 5 10^A/~^s~^, A — .Is^^, to — -Is, Rito — 5 10~''M. 
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